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ABSTRACT 

We present multi-epoch KPNO Coude Feed observations of interstellar Na I 
and Ca II absorption toward the bright star p Leo. Comparisons of the Na I 
profiles observed over a period of 8 years reveal significant temporal variations 
in the Na I column in at least one component, implying that there is "structure" 
at scales of order of the proper motion (~12 AU). Archival HST Goddard High 
Resolution Spectrograph observations of the C I fine-structure excitation in the 
variable component suggest that the density is <20 cm -3 , significantly lower than 
the densities inferred in past H I 21 cm and Na I studies. We suggest that the 
bulk of the trace neutral species are in the density peaks within an interstellar 
cloud. The patchy distribution of these species naturally gives rise to the large 
fluctuations seen on scales of 10-1,000 AU in past temporal and binary studies. 
This picture predicts that the scales over which fluctuations will be observed vary 
as a function of the ionization rate of a species. 

Subject headings: ISM: clouds - ISM: structure - stars: individual (HD 91316) 



1. Introduction 

The sightline toward p Leo (HD 91316) has long been a favorite of spectroscopists, 
even before the interstellar origin of the so-called "stationary" lines was established. Harper 



1 Based on observations with the NASA/ESA Hubble Space Telescope, obtained from the data archive at 
the Space Telescope Science Institute. STScI is operated by the Association of Universities for Research in 
Astronomy, Inc. under NASA contract NAS 5-26555. 

2 Visiting Astronomer, Kitt Peak National Observatory (KPNO). KPNO is a part of the National Optical 
Astronomy Observatories (NOAO), which is operated by the Association of Universities for Research in 
Astronomy , Inc. (AURA) under contract to the National Science Foundation (NSF). 
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(1914) first noticed that the Ca II K line toward p Leo showed the characteristic difference 
in velocity attributable to "stationary" lines, although the line was much wider than was 
typically observed with hints that there were two blended lines present. Beals (1936) showed 
that the Ca II lines toward p Leo were indeed complex, making it one of the first three 
sightlines (along with e and £ Ori) where multiple components were known. Subsequently, 
Sanford, Merrill, & Wilson (1938) showed that the Na I lines also have multiple components, 
while Adams (1943) showed that the "violet" Ca II component of Beals (1936) was itself 
multiple. The velocities and equivalent widths for the Na I and Ca II lines measured from 
photographic spectra by Routly & Spitzer (1952) are (to within the errors) the same as 
those measured today. Optical observations of p Leo have continued to the present day with 
a multitude of different instruments, and have shown the complex nature of the Na I, K I, 
and Ca II profiles (Hobbs 1969, 1971, 1974; Welty, Hobbs, & Kulkarni 1994; Welty, Morton, 
& Hobbs 1996; Welty & Hobbs 2001). In addition, observations of numerous ultraviolet 
absorption lines have been made by Copernicus, IUE, and i^ST/Goddard High Resolution 
Spectrograph (GHRS). 

Over the past decade, high spatial resolution optical and radio studies of the diffuse ISM 
have found strong evidence of pervasive subparsec-scale variations in interstellar absorption 
lines (Lauroesch 2001; Crawford 2003), although the precise interpretation of these fluctu- 
ations is not yet well understood (Heiles 1997; Elmegreen 1999; Deshpande 2000; Faison 
& Goss 2001; Welty & Fitzpatrick 2001; Lauroesch 2001; Crawford 2003). Initially, these 
observations were interpreted as ubiquitous concentrations of apparently dense atomic gas 
(n H > 10 3 cm -3 ) in otherwise diffuse sight-lines. However, recent direct measurements of 
the densities in a few of the optically selected clouds have suggested that in general the 
densities are much lower than previously inferred (Lauroesch et al. 1998; Lauroesch, Meyer, 
& Blades 2000; Welty & Fitzpatrick 2001), although observations in clouds identified by 
tracers of dense gas (K I) suggest that at least some of these features are associated with 
denser structures (Pan, Federman, & Welty 2001; Crawford 2002). Thus, it has been sug- 
gested that small-scale variations in the physical parameters may be responsible for much 
(if not all) of the observed variation in optical studies of small-scale structure. Recently, as 
part of a larger survey to identify very small scale (of order 10 AU) variations in interstellar 
lines, evidence for temporal (proper motion induced) variations in the interstellar Na I line 
toward p Leo have been observed. In this Letter, we will discuss these new observations 
as well as our analysis of archival Goddard High Resolution Spectrograph (GHRS) echelle 
observations of this sightline, and then discuss a simple model of neutral formation in the 
context of these results. 
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2. Observations and Data Reduction 

Rho Leo is a bright (V=3.8) Bllb runaway star, located at a distance of about 870 
pc approximately 690 pc above the galactic plane (Diplas & Savage 1994). The proper 
motion is 6.6±0.9 mas/yr (Perryman et al. 1997), which corresponds to a projected motion 
of 1.56±0.21 AU/yr. Since 1989, we have obtained multiple observations of the interstellar 
Na I and Ca II absorption lines toward p Leo using the KPNO Coude Feed and spectrograph 
at resolutions of ~1.4 and 3 km s _1 (see Table 1). These datasets have been reduced using 
the NOAO IRAF 3 data reduction packages in the same fashion as that used by Watson & 
Meyer (1996) and Lauroesch & Meyer (1999). The extracted continuum-fitted spectra for 
all observing runs are shown in Figure 1 - note in particular the variation in the Na I profile 
at a heliocentric velocity of 18 km s -1 . 

Table 1: p Leo Observation Log 



Species 


Date 


Resolution 


S/N ratio 






(km s _1 ) 


(per pixel) 


Nal 


11/95 


1.4 


370 




11/01 


3.4 


200 




11/02 


3.4 


220 




02/03 


1.4 


250 


Ca II 


05/89 


4.8 


190 




01/92 


3.6 


400 




11/02 


3.6 


330 


Al IIP 


12/1992 


3.6 


35-55 


C P 6 


03/1996 


3.5 


70 


Zn IP 


12/1992 


3.4-3.6 


70-80 



a Data obtained from the HST archive, obtained 

as parts of programs 2251 and 5882 (PI Hobbs). 
b Observations of the AA 1328A multiplet. 



In addition to the ground-based data, we have taken advantage of existing archival 
GHRS spectra of p Leo. Observations using multiple settings of the ECH-A and ECH-B 
gratings taken in March 1996 and December 1992 (respectively) are available from HST 
observing programs 5882 and 2251 (PI Hobbs). These datasets were reduced using the stan- 



3 IRAF is distributed by NOAO, which are operated by AURA, Inc. under cooperative agreement with 
the NSF. 
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Fig. 1.— Observations of the Na I D 2 (5889A), Na I Di (5895A), and Ca II K (3933A) lines 
towards p Leo. The top panel shows our two high resolution (1.4 km s -1 ) observations, note 
in particular the variation in the Na I profile near 18 km s _1 , as well as the somewhat more 
subtle variation near -1 km s _1 . Also shown at the bottom of the plot is the 2003 profile 
divided by the 1995 profile (with an offset of -0.9 applied to the relative intensity for plotting 
purposes). The lower panel shows our three moderate resolution observations of the Ca II 
K line. 
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dard STSDAS routines (POFFSETS, DOPOFF, and SPECALIGN) for both the WSCAN 
(ECH-B) and FP-SPLIT (ECH-A) datasets, and then selected transitions were continuum 
fitted and analyzed (Figure 2). 

The column densities, b-values, and relative velocities of the various components in both 
the KPNO and GHRS datasets were derived by profile fitting the continuum fitted spectra 
using the program f its6p (Welty, Hobbs, & York 1991). The wavelengths and oscillator 
strengths were taken from Morton (1991) and Welty et al. (1999), with the hyperfine 
splitting of the Na I D lines taken into account (Welty, Hobbs, & Kulkarni 1994). The 
Na I component model of Welty, Hobbs, & Kulkarni (1994) was used as a starting point for 
fitting the neutral species, with additional input from the Ca II component model of Welty, 
Morton, & Hobbs (1996) used for fitting the ionized species. The column density, b-value, 
and velocity of the various components were allowed to be free parameters in the initial fits 
- in the final fits the relative velocities of the components were fixed and the variation in 
the b-value limited for most components. 



3. The v~18 km s 1 Cloud 

As noted above, the most striking feature when comparing the Na I profiles in Figure 1 
is the strong temporal variation near 18 km s" 1 . Also significant, however, is the lack of 
a corresponding variation in the observed Ca II profiles, as well as the lack of convincing 
Na I profile variations at other velocities. We should note here that part of the reason for 
including p Leo in our survey was the apparent lack of H I 21-cm emission associated with 
the 18 km s" 1 component in the survey of Hobbs (1971), although more recent surveys show 
some 21-cm emission at this velocity in the region around p Leo (Hartmann & Burton 1997). 
In addition to showing clear evidence for temporal variability, the v~18 km s^ 1 component 
is also relatively isolated and so we will concentrate on this component for the rest of the 
discussion. 

The column densities for selected species derived for the v~18 km s _1 component are 
given in Table 2. The first thing to note is the large Na I differences between the years 
1989 and 2003, corresponding to a decrease of almost one-third in the column. However, 
there are two important caveats - first, we assumed the same b-value (0.58 km s -1 ) for the 
component as that found by Welty, Hobbs, & Kulkarni (1994) using their higher resolution 
data. Second, and perhaps more important, the published Welty, Hobbs, & Kulkarni (1994) 
profile shows no sign of the weak absorption redward of this component seen in our datasets. 
We note that this shelf of absorption is seen in the profile published by Hobbs (1971, 1974) 
which has a similar resolution to our datasets. This suggests that these features may have 
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Fig. 2.— The C I 1328 multiplet and the Zn II 2026 line toward p Leo, with the 2003 Na I 
observation shown for scale. Note the lack of excited C I* and C I** absorption - the plots 
are centered on the strongest transitions at 1329.1004A and 1329.5775A respectively. 
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been over-resolved in the Welty, Hobbs, & Kulkarni (1994) data and "lost" during continuum 
placement, perhaps altering the derived column density. In any case, as shown in Figure 1, 
there are significant changes in the Na I column even over the past 8 years, with a decrease 
of ~20% if we assume the Welty, Hobbs, & Kulkarni (1994) b-value. A period of 8 years 
corresponds to a projected proper motion at the distance to the star of ~ 12 AU. One useful 
comparison that can be made immediately is to compare the N(H) one would estimate from 
the observed Na I columns, to that which would be derived using the observed N(Zn II). The 
average of the 1989 and 1995 Na I columns is 27.5 xlO 10 cm -2 , from which one would infer an 
N(H) of ~5.5xl0 19 cm' 2 (Ferlet, Vidal-Madjar, & Gry 1995). Similarly, one would estimate 
N(H)~4.6-7.2xl0 19 cm" 2 from N(Zn II) (assuming a Zn depletion of 0.2 to 0.4 dex (Welty 
et al. 1999)). This consistency immediately highlights a potentially important difference 
between this cloud and the clouds studied previously in the ultraviolet by Lauroesch et al. 
(1998) and Lauroesch, Meyer, & Blades (2000), namely the observed N(Zn II) (and inferred 
N(H)) is approximately an order of magnitude larger for this cloud, even though the N(Na I) 
is only a factor of 2-3 larger. 

Table 2: Selected Column Densities at v=18 km s _1 toward p Leo 



Species 


Year 


Lines/Mult iplet 


N(X) 






Used (A) 


(cm" 2 ) 


Nal 


1989 a 


5895 


30.7xl0 10 




1995 


5889,5895 


24.3±0.5xl0 10 b 




2001 


5889,5895 


20.8±1.0xl0 10 b 




2002 


5889,5895 


20.2±1.0xl0 10 b 




2003 


5889,5895 


20.9±0.6xl0 10 b 


Al III 


1992 


1854, 1862 


5.9±l.lxlO n 


C I 


1996 


1328 


57.7±6.2xlO n 


C I* 


1996 


1328 


<8.4xlO n d 


c i** 


1996 


1328 


<8.6xlO n d 


Cr II 


1992 


2056, 2062 


1.8±0.2xlO n 


Mg I 


1992 


2026 


17.6±1.8xlO n 


Zn II 


1992 


2026, 2062 


13.1±1.8xlO n 



We can now use the other species listed in Table 2 to estimate the physical conditions in 
this cloud. We can use the observed N(C I) and the upper limit on N(C I*) to place limits on 
the pressure and density in this cloud (Jenkins & Shaya 1979) - our upper limits to N(C I*) 
imply for T~100-500 K (consistent with the observed b-value) pressures of log(p/k)<3.3-3.8 
cm -3 K. For an ideal gas at these temperatures, such pressures correspond to densities of 
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n(H)<20 cm -3 , similar to the densities found in previous studies (Lauroesch et al. 1998; 
Lauroesch, Meyer, & Blades 2000). In addition, the fine-structure equilibrium can also be 
used to estimate the minimum distance between the cloud and p Leo itself, since photon- 
pumping can significantly alter the observed level populations (Jenkins & Shaya 1979; 
Lauroesch et al. 1998). The observed upper limits to N(C I*) suggest that this cloud is >10 
pc away from p Leo. The upper limit on the density coupled with our estimates of N(H) 
above suggest a cloud size of order 1 pc, much larger than the length scale over which we 
are seeing Na I fluctuations (~12 AU) but also smaller than the estimated distance between 
the cloud and p Leo. The observed ratio of N(Cr II)/N(Zn II)~1.4 would suggest that this 
cloud has a "warm cloud" depletion pattern (Welty et al. 1999). If we now assume a typical 
"warm cloud" value for the depletion of Al, we would estimate N(A1 II+A1 III)~l.lxl0 13 
cm -2 from N(Zn II), which would suggest that while Al II is the dominant ionization state 
in this cloud there is likely to be some degree of photo-ionization. We can also use the 
measured neutral and Zn II column densities to estimate electron densities of 0.04, 0.15, and 
0.05 cm -3 from C I, Na I and Mg I respectively assuming T=100 K and a 'warm cloud" 
depletion pattern (Pequignot & Aldrovandi 1986; Lauroesch, Meyer, & Blades 2000). If we 
assume the dominant source of electrons is singly ionized carbon, these values correspond to 
a density of ~100-400 cur 3 at T=100 K, significantly in excess of the upper limit derived 
above. 



4. Biased Neutral Formation in Interstellar Gas? 



The parameters derived above for the v~18 km s" 1 component toward p Leo appear 
fairly typical for a diffuse cloud. The column density of hydrogen is ~6xl0 19 cm" 2 , and 
we would estimate a cloud size of ~1 pc along the line of sight. The density of hydrogen 
appears to be of order 10 cm -3 , with the derived electron density and observed Al III 
suggesting some partial ionization of hydrogen. There is one significant problem that still 
remains however, namely the origin of the observed fluctuations in the Na I profile over 
relatively short timescales. Since the photo-ionization and recombination rates for Na I 
are relatively low (Pequignot & Aldrovandi 1986), one would not expect to see significant 
variations in the column on timescales shorter than thousands of years in an interstellar 
cloud. If we were thus to assume the observed fluctuations set the transverse size of the 
cloud, the earliest observations of p Leo suggest that significant Na I and Ca II absorption 
has been detectable at this velocity for at least 65 years, setting a minimum transverse size of 
100 AU for this "structure" (at the distance of the star). The result is similar to the models 
suggested by Heiles (1997), namely a very long, thin sheet or filament with a thickness 
of order one thousandth the length. However, there is an alternative explanation for this 
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observed variation. 

Consider an isolated diffuse interstellar cloud - it will have within it density fluctuations 
at scales much smaller than a parsec (Elmegreen 1999; Cho, Lazarian, & Vishniac 2002). 
Just as the force of gravity biases galaxies to form in the peaks of the cosmic density distribu- 
tion, the neutral species in interstellar clouds are biased by the recombination rates to form 
on the peaks of the density distribution within a cloud. It is well known that recombination 
rates have a strong density dependence (Pequignot & Aldrovandi 1986), indeed this formed 
the basis for previous suggestions that relatively small density fluctuations could give rise to 
the observed small scale structure (Lauroesch & Meyer 1999; Lauroesch, Meyer, & Blades 
2000). Furthermore, the C I survey of Jenkins & Tripp (2001) has identified a population of 
high pressure components in the ISM, which could represent these density peaks. We now 
ask what are some of the observational consequences of such a picture? 

First we can estimate a minimum apparent length scale for the size of fluctuations 
observable in any trace species. The scale is just set by the distance a neutral can travel 
before it can be ionized, which in low density gas is of order the velocity of an atom divided 
by the photo-ionization rate. Assuming v~0.5 km s~ x which is a "typical" b-value for the 
known temporally variable components, we derive an apparent scale of order 250 AU for 
Na I structures. Despite our having neglected to correct for collisions, this size is of order 
that observed in temporal fluctuations to date, and is consistent with the observations of 
Na I towards binary star systems. Furthermore, note that this scale is dependent upon the 
photo-ionization rate for a given species, and thus (for example) for Ca II we would obtain 
a scale about 5 times larger. This result very simply explains why we do not observe a 
corresponding variation in the Ca II profile towards p Leo. In addition, it suggests that 
species with higher photo-ionization rates such as Ca I and (perhaps) C I should show larger 
column density variations than species such as Na I or K I. Indeed, Crawford (2002) used 
observations of Ca I to suggest that the variable component towards k Velorum arises in 
very dense material (nn > 10 3 cm~ 3 ). 

A further consequence of this picture is that only a fraction of the total gas would reside 
in these peaks. Thus, adjacent sightlines should show smaller changes in the columns of the 
dominant ions, since the size of the peaks is only a fraction of the size of a cloud the additional 
column density for the dominant species would be significantly reduced. Consequently, the 
neutral fraction within the density peaks would be higher (perhaps much higher) than what 
one would calculate by taking the integrated columns, suggesting that the electron densities 
inferred from ionization balances are too low. Furthermore, the fine-structure equilibrium of 
C II and Si II would also be affected. Just like the neutral species, most of the observed C II* 
and Si II* absorption would arise in the density peaks, but the columns of the ground-states 
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would reflect the large scale distribution of gas and would include the perhaps dominant 
contribution by the material outside the peaks. Thus, in this picture, these methods provide 
only an ill-defined "average" electron density within a cloud, and would explain why there 
is so much difficulty in deriving self-consistent electron densities in many cases (Welty et al. 
1999). Assumptions about partial ionization of hydrogen could then not be made using the 
(separately) derived electron and total hydrogen densities, since these tracers will likely be 
probing different portions of a cloud. 

A final question is whether there is a connection between the observed small scale 
variations in Na I and the high pressure clumps observed by Jenkins & Tripp (2001). In this 
picture, the fine-structure equilibrium of C I becomes a measure of some sort of "average" 
density of the peaks along a sightline through a cloud. To date, there have been limited 
observations of the C I fine structure equilibrium associated with optically selected small scale 
structures (Lauroesch et al. 1998; Lauroesch, Meyer, & Blades 2000; Welty & Fitzpatrick 
2001). Including p Leo, the observed Na I structures in all of these sightlines appear to 
arise in relatively low density gas, and not in the high pressure clumps detected by Jenkins 
& Tripp (2001). However, as Crawford (2003) recently noted for the temporally variable 
components towards both HD 32040 and HD 219188, these sightlines may be sampling a 
mixture of low pressure gas with a fraction (perhaps 20%) of high pressure, dense material 
with which the temporal variation may be associated. In the picture developed above, the 
observed temporal and binary star variations are arising in the neutral "halo" around one 
or more of these high pressure clumps, implying that these clumps are a fraction of the 
apparent size we calculated above. Yet another possibility is that we are seeing the remnant 
neutrals formed in transient high pressure clumps which have already expanded into the 
ambient cloud. Planned re-observations of both HD 32040 and HD 219188 will likely help in 
beginning to resolve this question. In any case, this suggests that the density of the majority 
of the gas in a cloud is smaller than that inferred from the C I observations, but that some 
fraction resides in higher density structures within which arises the bulk of the trace neutral 
species. 

The authors thank the referee Ian Crawford for his valuable comments on this manuscript. 
It is also a pleasure to acknowledge the support of the staff of KPNO, especially Daryl Will- 
marth, for their assistance in obtaining this data. 
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